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Abstract. The semiclassical limit of a weakly coupled nonlinear focusing Schrodinger sys- 
^\ . tem in presence of a nonconstant potential is studied. The initial data is of the form (wi, U2) 

I with Ui = ri('^^)ei^'^, where (ri,r2) is a real ground state solution, belonging to a suit- 

able class, of an associated autonomous elliptic system. For e sufficiently small, the solution 
(01, (j)2) will been shown to have, locally in time, the form (ri |^ ^~^(*) ^g7^-g(*)^ ?'2( ^~^*-*'' )e'^'^^*-'), 
where {x{t),£,{t)) is the solution of the Hamiltonian system x{t) ~ ^(t), £^{t) ~ ~'S/V{x{t)) 
^ ■ with x{0) = i and ^(0) = |. 

O 

(N 

1. Introduction and main result 

' 1.1. Introduction. In recent years much interest has been devoted to the study of systems 

■ of weakly coupled nonlinear Schrodinger equations. This interest is motivated by many 

J3 ' physical experiments especially in nonlinear optics and in the theory of Bose-Einstein con- 

densates (see e.g. [1, 17, 24, 26]). Existence results of ground and hound states solutions have 
been obtained by different authors (see e.g. [3, 5, 13, 21, 22, 30]). A very interesting aspect 
^ . regards the dynamics, in the semiclassical limit, of a general solution, that is to consider 

. the nonlinear Schrodinger system 
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^e^t(j)l + ^A0i - V{x)<j), + Mlhl^" + f3\<l>2r'\hr') = in x 

tedt(j)2 + yA02 - V{x)^2 + 02(|02|'" + /3|0ir+V2r~') = iu X 
0i(O,x) = </.?(x) M0,x) = 4{x), 



^ . with 0<p<2/iV, A^>1 and /5 > is a constant modeling the birefringence effect of the 

^ . material. The potential V{x) is a regular function in modeling the action of external 

forces (see (1.11)), (pi : x C are complex valued functions and e > is a small 

parameter playing the role of Planck's constant. The task to be tackled with respect to this 
system is to recover the full dynamics of a solution (</)f , </)|) as a point particle subjected to 
galileian motion for the parameter e sufficiently small. Since the famous papers [2, 14, 16], 
a large amount of work has been dedicated to this study in the case of a single Schrodinger 
equation and for a special class of solutions, namely standing wave solutions (see [4] and the 



2000 Mathematics Subject Classification. 34B18, 34G20, 35Q55. 

Key words and phrases. Weakly coupled nonlinear Schrodinger systems, concentration phenomena, semi- 
classical limit, orbital stability of ground states, soliton dynamics. 

The first and the second author are supported by the MIUR national research project "Variational 
Methods and Nonlinear Differential Equations" , while the third author is supported by the 2007 MIUR 
national research project "Variational and Topological Methods in the Study of Nonlinear Phenomena" . 

1 



2 



E. MONTEFUSCO, B. PELLACCL AND M. SQUASSINA 



references therein). When considering this particular kind of solutions one is naturally lead 
to study the following elliptic system corresponding to the physically relevant case p = 1 
(that is Kerr nonlinearities) 



(1.2) 



-e'^Au + V{x)u = u^ + (3v'^u in M^, 
-e'^Av + V{x)v = v^ + (3u^v in M^, 



so that the analysis reduces to the study of the asymptotic behavior of solutions of an elliptic 
system. The concentration of a least energy solution around the local minima (possibly 
degenerate) of the potential V has been studied in [27], where some sufficient and necessary 
conditions have been established. To our knowledge the semiclassical dynamics of different 
kinds of solutions of a single Schrodinger equation has been tackled in the series of papers [7, 
18, 19] (see also [6] for recent developments on the long term soliton dynamics), assuming 
that the initial datum is of the form r((x — x)/e)e^^'^, where r is the unique ground state 
solution of an associated elliptic problem (see equation (1.8)) and ^ e M^. This choice of 
initial data corresponds to the study of a different situation from the previous one. Indeed, 
it is taken into consideration the semiclassical dynamics of ground state solutions of the 
autonomous elliptic equation once the action of external forces occurs. In these papers it is 
proved that the solution is approximated by the ground state r-up to translations and phase 
changes-and the translations and phase changes are precisely related with the solution of 
a Newtonian system in governed by the gradient of the potential V. Here we want to 
recover similar results for system (1.1) taking as initial data 



;i.3) 0?(x) = nf^)eM, ^l(x) = rJ^)ei-« 



where the vector R = (ri, r2) is a suitable ground state (see Definition 1.3) of the associated 
elliptic system 

1 



--An + n = n{\n\'p + /3|r2|p+>i|p-') in 

--Ar2 + r2 = radral^P + P\ri\P+^\r2\''~^) in 



2 

When studying the dynamics of systems some new difficulties can arise. First of all, we 
have to take into account that, up to now, it is still not known if a uniqueness result (up 
to translations in M^) for real ground state solutions of (E) holds. This is expected, at 
least in the case where j3 > 1. Besides, also nondegeneracy properties (in the sense provided 
in [12, 28]) are proved in some particular cases [12, 28]. These obstacles lead us to restrict the 
set of admissible ground state solutions we will take into consideration (see Definition 1.3) 
in the study of soliton dynamics. 

Our first main result (Theorem 1.5) will give the desired asymptotic behaviour. Indeed, 
we will show that a solution which starts from (1.3) (for a suitable ground state R) will 
remain close to the set of ground state solutions, up to translations and phase rotations. 
Furthermore, in the second result (Theorem 1.9), we will prove that the mass densities 
associated with the solution (pi converge-in the dual space of C'^{W'^) x C^(]R^)-to the delta 
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measure with mass given by ||ri||i2 and concentrated along x{t), solution to the (driving) 
Newtonian differential equation 

(1.4) x{t) = -Vy(x(t)), x{0) = X, x{0) = i 

where x and ^ are fixed in the initial data of (1.1). A similar result for each single component 
of the momentum density is lost as a consequence of the birefringence effect. However, we 
can afford the desired result for a balance on the total momentum density. This shows that- 
in the semiclassical regime-the solution moves as a point particle under the galileian law 
given by the Hamiltonian system (1.4). In the case of V constant our statements are related 
with the results obtained, by linearization procedure, in [31] for the single equation. Here, 
by a different approach, we show that (1.4) gives a modulation equation for the solution 
generated by the initial data (1.3). Although we cannot predict the shape of the solution, 
we know that the dynamic of the mass center is described by (1.4). The arguments will 
follow [7, 18, 19], where the case of a single Schrodinger equations has been considered. The 
main ingredients are the conservation laws of (1.1) and of the Hamiltonian associated with 
the ODE in (1.4) and a modulational stability property for a suitable class of ground state 
solutions for the associated autonomous elliptic system [E), recently proved in [28] by the 
authors in the same spirit of the works [31, 32] on scalar Schrodinger equations. 

The problem for the single equation has been also studied using the WKB analysis (see 
for example [9] and the references therein), to our knowledge, there are no results for the 
system using this approach. Some of the arguments and estimates in the paper are strongly 
based upon those of [19]. On the other hand, for the sake of self-containedness, we prefer 
to include all the details in the proofs. 

1.2. Admissible ground state solutions. Let be the space of the vectors $ = (0i, 02) 
in H = if^(M^; C^) endowed with the rescaled norm 

\mk = ^mi + ^j'^ni 

where = ||(0i,02)||2 = ll^illi + ll'/'slli and ||0i||2 = 110^11^2 is the standard norm in the 

Lebesgue space given by ||0i||2 = / 0i(a;)0i(x)(ix. 

We aim to study the semiclassical dynamics of a least energy solution of problem (E) once 
the action of external forces is taken into consideration. 

In [3, 22, 30] it is proved that there exists a least action solution R = (ri,r2) ^ (0,0) 
of (E) which has nonnegative components. Moreover, is a solution to the following 
minimization problem (cf. [23, Theorems 3.4 and 3.6]) 

(1.5) S{R)=mmS, where M = {U e M : llUh = WRh} , 

M 

where the functional £^ : EI ^ M is defined by 

(1.6) S{U) = ^\\VU\\l- I F^{U)dx 

(1.7) F,{U) = |n2p^+2 + 2/3|nir>2r^), 
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for any U = (^1,^2) G H. We shall denote with Q the set of the (complex) ground state 
solutions. 

Remark 1.1. Any element V = {vi, V2) of Q has the form 

V{x) = {e'^'\vi{x)\,e'^'\v2{x)\), x E M^, 

for some 9i,92 G (so that (|t'i|, \v2\) is a real, positive, ground state solution). Indeed, if 
we consider the minimization problems 

ac = inf {S{V) -.VeU, \\V\\l2 = \\R\\l^}, 

aK = mi{S{V) : V e H\R^;R^) \\V\\l2 = WRWl^} 

it results that ac = ctr. Trivially one has ac < ctr- Moreover, ii V = (fi,f2) G H, due to 
the well-known pointwise inequality |V|wi(x)|| < |Vfi(a;)| for a.e. x G M^, it holds 



so that also i^dt'i], \v2\) < ^iV). In particular, we conclude that cr^ < ac, yielding the de- 
sired equality ac = cr^- Let now V = (f 1, f 2) be a solution to ac and assume by contradiction 
that, for some i = 1,2, 

C^{{x G : \V\vi\{x)\ < \Vvi{x)\}) > 0, 

where is the Lebesgue measure in M^. Then ||(|t>i|, 11121) ||l2 = ||^||l2, and 
2 2 
cr]8<^^y" \V\vi\\^dx- j Fp{\vi\,\v2\)dx <]^'Y^ j \Vvi\^dx- j Fp{vi,V2)dx = ac, 

which is a contradiction, being ac = ajg. Hence, we have |V|wj(x)|| = |Vwj(x)| for a.e. 
X G and any i = 1,2. This is true if and only if TZeViViXmvi) = TmviV (JZevi). In 
turn, if this last condition holds, we get 

ViWvi = IZeViW [TZeVi] + XmViV{ImVi), a.e. in M^, 

which implies that TZe {wi{x)'Vvi{x)) = a.e. in R^. Finally, for any i = 1,2, from this last 
identity one immediately finds 6i G with Vi = e*^'|wj|, concluding the proof. 

In the scalar case, the ground state solution for the equation 

(1.8) _ Ia^ + r- = r^P+i inM^ 

^ ' 2 

is always unique (up to translations) and nondegenerate (see e.g. [20, 25, 31]). For sys- 
tem (E), in general, the uniqueness and nondegeneracy of ground state solutions is a delicate 
open question. 

The so called modulational stability property of ground states solutions plays an important 
role in soliton dynamics on finite time intervals. More precisely, in the scalar case, some 
delicate spectral estimates for the seld-adjoint operator £"{r) were obtained in [31, 32], 
allowing to get the following energy convexity result. 
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Theorem 1.2. Le r be a ground state solution of equation (1.8) with p < 2/N. Let (f) G 

iJ^(]R^,C) be such that \\4>\\2 = \\r\\2 (^nd define the positive number 

r^= inf l|0(-)-e^^r(--y))||i.. 
Then there exist two positive constants A and C such that 

r^<c(£(0)-^(i?)), 

provided that S{(f>) — £{R) < A. 

For systems, we consider the following definition. 

Definition 1.3. We say that a ground state solution R = (ri, r2) of system (E) is admissible 
for the modulational stability property to hold, and we shall write that R ifri G if^(M^) 
are radial, \x\ri G L^(R^), the corresponding solution 4>i{t) belongs to iir^(R^) for all times 
t > and the following property holds: let ^ eM be such that ||$||2 = ll-^lh and define the 
positive number 

(1.9) U:= inf l|$(-)-(e^V(--l/),e^'V2(--2/))||^. 

6li,6»2e[0,27r) 

Then there exist a continuous function p : Ms^ Ms^ with as ^ 0^ and a positive 

constant C such that 

p(r$) + r$ < c{S{<^)-s{R)). 

In particular, there exist two positive constants A and C such that 

(1.10) <C'{S{<^) -S{R)), 
provided that < A. 

In the one dimensional case, for an important physical class, there exists a ground state 
solution of system (E) which belongs to the class TZ (see [28]). 

Theorem 1.4. Assume that N = 1, p E [1, 2) and (3 > 1. Then there exists a ground state 
solution R = (ri,r2) of system (E) which belongs to the class TZ. 

1.3. Statement of the main results. The action of external forces is represented by a 
potential V : M satisfying 

(1.11) V is a, function bounded with its derivatives, 

and we will study the asymptotic behavior (locally in time) as e ^ of the solution of the 
following Cauchy problem 

.2 



tedth + yA0i - V{x)(P, + 0i(|0ip^ + /^l^ar+ViT"') = in x R+, 
^e^t(f)2 + -A02 - ^(x)02 + 02(|02|'^ + /3|0ir+V2r"') = in X M+, 
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where x, ^ G N > 1, the exponent p is such that 
(1.12) 0<p<2/A^ 

It is known (see [15]) that, under these assumptions, and for any initial datum in L^, there 
exists a unique solution $^ = {(j)'i,(j)2) of the Cauchy problem that exists globally in time. 
We have chosen as initial data a scaling of a real vector R = (ri, r2) belonging to TZ. 

The first main result is the following 

Theorem 1.5. Let R = (ri,r2) be a ground state solution of (E) which belongs to the class 
TZ. Under assumptions (1.11), (1-12), let = (0^,02) the family of solutions to system 
(^e). Furthermore, let {x{t),^{t)) be the solution of the Hamiltonian system 



;i.i3) 



x{t) = m 
m = -vv(x(t)) 

x(0) = X 

{m = l 



Then, there exists a locally uniformly bounded family of functions 9f : S^, i = 1,2, 

such that, defining the vector Q^it) = (gf (x, t), g|(x, t)) by 



it holds 

(1.14) \\^%t) - Qe{t)U < 0{e), ase^O 

locally uniformly in time. 

Roughly speaking, the theorem states that, in the semiclassical regime, the modulus of 
the solution is approximated, locally uniformly in time, by the admissible real ground 
state (ri,r2) concentrated in x{t), up to a suitable phase rotation. Theorem 1.5 can also be 
read as a description of the slow dynamic of the system close to the invariant manifold of 
the standing waves generated by ground state solutions. This topic has been studied, for 
the single equation, in [29]. 

Remark 1.6. Suppose that ^ = and x is a critical point of the potential V. Then the 
constant function {x{t),C,{t)) = {x, 0), for all t G M"*", is the solution to system (1.13). As a 
consequence, from Theorem 1.5, the approximated solutions is of the form 



\ e / 

that is, in the semiclassical regime, the solution concentrates around the critical points of 
the potential V. This is a remark related to [27] where we have considered as initial data 
ground states solutions of an associated nonautonomous elliptic problem. 
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Remark 1.7. As a corollary of Theorem 1.5 we point out that, in the particular case of a 
constant potential, the approximated solution has components 

'X — X — ^t" 



Hence, the mass center x{t) of moves with constant velocity ^ realizing a uniform 

motion. This topic has been tackled, for the single equation, in [31]. 

Remark 1.8. For values oi (3 > 1 both components of the ground states R are nontrivial 
and, for R & TZ, the solution of the Cauchy problem are approximated by a vector with 
both nontrivial components. We expect that ground state solutions for /3 > 1 are unique 
(up to translations in M^) and nondegenerate. 

We can also analyze the behavior of total momentum density defined by 

(1.15) P%x,t) ■=pl{x,t) +pl{x,t), forxGM^, t>0, 
where 

(1.16) pt{x,t) := ^Jm(0^(x,t)V0f(x,t)), for z = 1, 2, x G M^, t > 0. 

Moreover, let M{t) := (mi + m2)^{t) be the total momentum of the particle x{t) solution 
of (1.13), where 

(1.17) mj := ||rj||2, fori = 1,2. 

The information about the asymptotic behavior of and of the mass densities l^fp/e^ 
are contained in the following result. 

Theorem 1.9. Under the assumptions of Theorem 1.5, there exists Eq > such that 
||(|0^|V5^rfx, m'/e''dx) -(mi,m2)5.(i)||(^,^^,^, < 0{e'), 
\\P%t,x)dx-M{t)6,(^t)\\^^,^, < 0{e^), 
for every e G {0,eo) and locally uniformly in time. 

Remark 1.10. Essentially, the theorem states that, in the semiclassical regime, the mass 
densities of the components (pi of the solution $^ behave as a point particle located in x{t) of 
mass respectively nii and the total momentum behaves like M{t)Sx{t)- It should be stressed 
that we can obtain the asymptotic behavior for each single mass density, while we can only 
afford the same result for the total momentum. The result will follow by a more general 
technical statement (Theorem 2.4). 

Remark 1.11. The hypotheses on the potential V can be slightly weakened. Indeed, we 
can assume that V is bounded from below and that d°'V are bounded only for |a| = 2 or 
|a| = 3. This allows to include the important class of harmonic potentials (used e.g. in 
Bose- Einstein theory), such as 

1 ^ 
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Hence, equation (1-13) reduces to the system of harmonic oscillators 
(1.18) x,{t)+uj]x,{t) = 0, j = l,...,iV. 

For instance, in the 2D case, renaming xi(t) = x{t) and X2{t) = y{t) the ground states 
solutions are driven around (and concentrating) along the lines of a Lissajous curves having 
periodic or quasi-periodic behavior depending on the case when the ratio uo-Juj is, respec- 
tively, a rational or an irrational number. See Figures 1 and 2 below for the corresponding 
phase portrait in some 2D cases, depending on the values of Ui/uj. 




-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 



x(l) x(t) 



Figure 1. Phase portrait of system (1.18) in 2D with uji/uj2 = 3/5 (left) and 
ujxluji = 7/5 (right). Notice the periodic behaviour. 




Figure 2. Phase portrait of system (1.18) in 2D with uji/u2 = v^/3 increas- 
ing the integration time from t G [0,407r] (left) to t G [0,607r] (right). Notice 
the quasi-periodic behaviour, the plane is filling up. 



The paper is organized as follows. 
In Section 2 we set up the main ingredients for the proofs as well as state two technical 
approximation results (Theorems 2.2, 2.4) in a general framework. In Section 3 we will 
collect some preliminary technical facts that will be useful to prove the results. In Section 4 
we will include the core computations regarding energy and momentum estimates in the 
semiclassical regime. Finally, in Section 5, the main results (Theorems 1.5 and 1.9) will be 
proved. 
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2. A MORE GENERAL SCHRODINGER SYSTEM 

In the following sections we will study the behavior, for sufficiently small e, of a solution 



h,4>2) of the more general Schrodinger system 



iedt(f)2 + -7r^<P: 



0i(O,x) 



ri 



2 -iy(x)02 + 02(1021'^ + 



|P+1| 



|P+1| 



in 
= in 



pN 



pN 



X 



X 



e 



fx — X 

02(0, = r2( 



where p verifies (1.12), the potentials V, W both satisfy (1.11) and (ri,r2) is a real ground 
state solution of problem (E). As for the case of a single potential, we get a unique globally 
defined $^ = (01,01) that depends continuously on the initial data (see, e.g. [15, Theorem 
1]). Moreover, if the initial data are chosen in x H^, then $^(t) enjoys the same regularity 
property for all positive times t > (see e.g. [10]). 

Remark 2.1. With no loss of generality, we can assume V,W > 0. Indeed, if 0i,02 is a 
solution to (-Fe), since V, W are bounded from below by (1.11), there exist fi > such that 



he 



V{x) + /i > and W{x) + /i > 0, for all x G M^. Then 0i 
solution of (Fs) with V + ft (resp. W + ^) in place of V (resp. W). 



and 



926 



IS a 



We will show that the dynamics of 
X = (xi,X2) : 
of the following Hamiltonian systems 



^'!,0^2; 

T)27V 



Ut) = -VV{x^{t)) 
(xi(0),6(0)) = (£,fi) 



is governed by the solutions 
S = (ei,6):K^M27V, 



Ut) = -VW{x2{t)) 
(X2(0),6(0)) = (£,6) 



Notice that the Hamiltonians related to these systems are 



(2.1) 



and are conserved in time. Under assumptions (1.11) it is immediate to check that the 
Hamiltonian systems {H) have global solutions. With respect to the asymptotic behavior 
of the solution of (F^) we can prove the following results. 

2.1. Two more general results. We now state two technical theorems that will yield, as 
a corollary. Theorems 1.5 and 1.9. 



!,0!) he 



Theorem 2.2. Assume (1.12) and that V, W both satisfy (l.U). Let $^ = 

the family of solutions to system (-Fe). Then, there exist Eq > 0, > 0, a family of 
continuous functions g'^ : M"*" M with g'^{0) = 0{e'^), locally uniformly bounded sequences 
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of functions 6f : and a positive constant C , such that, defining the vector Qe{t) 

{qt{x,t), by 



qKx,t)=rI^^^^]e^^-^' 



e 

it results 



1,2 



Wit)-Q,it)\\u,<C^g^it) + 

for all e G (0,£:o) o.'^d all t G [0,T^], where Xi{t) is the first component of the Hamiltonian 
system for V in {H). 

Remark 2.3. Theorem 2.2 is quite instrumental in the context of our paper, as we cannot 
guarantee in the general case of different potentials that the function is small as e vanishes, 
locally uniformly in time. Moreover, the time dependent shifting of the components qi into 
Xi{t) is quite arbitrary, a similar statement could be written with the component X2{t) in 
place of this arbitrariness is a consequence of the same initial data x in {H) for both 

Xi and X2- The task of different initial data in [H) for xi and X2 is to our knowledge an 
open problem. 

In the following, if are the second components of the systems in (iJ), we set 

(2.2) M{t) := miei(t) + mg^W, t > 0. 

If = {(pl, 01) is the family of solutions to (-Fe), we have the following 

Theorem 2.4. There exist Eq > and > and a family of continuous functions g'^ : 
M+ ^ M with g%0) = 0{e^) such that 

\\P'{t,x)dx - M{t)6^,^t))\\^c^^, ^ 
for every e G (0, Sq) and all t G [0, T^]. 

3. Some preliminary results 

In this section we recall and show some results we will use in proving Theorems 1.5, 1.9, 
2.2 and 2.4. First we recall the following conservation laws. 



Proposition 3.1. The mass components of a solution $ of (F^), 

1 

e 

are conserved in time. Moreover, also the total energy defined by 



(3.1) Mnt) := ^mmh, /orz = l,2, t>0. 



(3.2) E%t) = El{t) + Em 
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is conserved as time varies, where 

Em 
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Proof. This is a standard fact. For the proof, see e.g. [15]. ■ 

Remark 3.2. From the preceding proposition we obtain that, due to the form of our initial 
data, the mass components M^it) do not actually depend on e. Indeed, for i = 1, 2, 



(3.3) 



\^dx = 



X — X 



dx = rrij. 



Thus, the quantities ^f/e^^^ have constant norm in equal, respectively, to m^. In The- 
orem 2.4 we will show that, for sufficiently small values of e, the mass densities behave, 
point-wise with respect to t, as a 5 functional concentrated in xi{t). 

In the following we will often make use of the following simple Lemma. 

Lemma 3.3. Let A G C^(M^) be such that A, DjA, DfjA are uniformly bounded and let 
R = (ri,r2) be a ground state solution of problem E. Then, for every y G fixed, there 
exists a positive constant Co such that 



(3.4) 



[A{ex + y) — A{y)] rf{x)dx 



< Cos'. 



Proof. By virtue of the regularity properties of the function A and Taylor expansion The- 
orem we get 



[A{ex + y) — A{y)] r1{x)dx 



<\\VA{y)\ 
+ ||Hes(A)||, 



xrf(x)dx 
x\'^r^{x)dx 



where ||Hes(y4)||oo denotes the L°° norm of the Hessian matrix associated to the function 
A. The first integral on the right hand side is zero since each component is radial. The 
second integral is finite, since \x\ri G L^(]R^). ■ 

In order to show the desired asymptotic behavior we will use the following property of 
the functional 6y on the space C^(]R^). 

Lemma 3.4. There exist Kq, Ki, K2 positive constants, such that, if \\5y — 5z\\c'^-'< Kq then 

Ki\y - z\ < \\6y - 6z\\c2* < K2\y - z\ 
Proof. For the proof see [19, Lemma 3.1, 3.2]. ■ 

The following lemma will be used in proving our main result. 
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Lemma 3.5. Let $^ = (0f,0|) be a solution of {F^) and consider the vector functions 
ai-.R^R^ defined by 

(3.5) ai{t) = J pl{x,t)dx -mi^i{t), t>0,i = l,2, 

where the ^iS are defined in {H) and the m^s are defined in (1.17), for i = 1,2. Then 
{t C(i{t)} is a continuous function and (0) = 0, for i = 1,2. 

Remark 3.6. The integral in (3.5) defines a vector whose components are the integral of 
Jm(0f90f/9xj)/e^-i for j = 1, . . . , N, so that : R ^ . 



Proof. The continuity of immediately follows from the regularity properties of the solu- 
tion In order to complete the proof, first note that, for all x G M^, 

SO that, real function, the conclusion follows by a change of variable. 



Lemma 3.7. Let V and W both satisfying assumptions (1.11) and let $^ = (</>i,0|) be a 
solution of (-Fe). Moreover, let A a positive constant defined by 

(3.6) A = K, sup [|xi(t)| + \x2it)\] + Ko 

[0,To] 

where Xi{t) is defined in {H), Kq and Ki are defined in Lemma 3.4, and let x be a C°°{R^) 
function such that < x ^ 1 and 

(3.7) x{x) = 1 if \x\ < A, x{x) = if \x\ > 2A. 
Then the functions 



rilit) = miV{xi{t)) - / xix)V{x)\(f)l{x,t)\^dx, 
rilit) = m2W{x2{t)) - / x{x)W{x)W2{x,t)\^dx. 



are continuous and satisfy |77|(0)| = 0{e'^) fori = 1,2. 

Proof. The continuity of r]f immediately follows from the regularity properties of the solu- 
tion We will prove the conclusion only for ?7i(0), the result for ?7|(0) can be showed in 
an analogous way. We have 



l^i(0)| 



rmV{x,{0))-^ I x{x)V{xM{x,0)\'dx 



miV{x) ^ / V{x)rly jdx 

+ ^ [ {l-x{x))V{x)rlf^)dx. 

^ J\x\>A \ e / 
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Then, by Lemma 3.3, and a change of variables imply 

|r7^(0)| < 0{e^) + j{l-x{x + ey)) V{x + ey)rl {y) dy. 

The properties of x ^"^^ ^ind assumption (1.11) yield the conclusion. ■ 

We will also use the following identities. 
Lemma 3.8. The following identities holds for i = 1,2. 

(3.9) —J^(x,t) = -dw,ptix,t), X G M^, t > 0. 

Moreover, for all t > 0, it results 



(3.10) / -^ix,t)dx = -^ / Vr(x)|0^(x,t)|2rfx-^ / \/Wixmix,t)fdx, 



where P'^{x,t) is the total momentum density defined in (1.15). 

Remark 3.9. It follows from identity (3.10) that for systems with constant potentials the 
total momentum J P^dx is a constant of motion. 

Remark 3.10. As evident from identity (3.10) as well as physically reasonable, in the case 
of systems of Schrodinger equations, the balance for the momentum needs to be stated for 
the sum P"^ instead on the single components . See also identities (3.11) and (3.12) in the 
proof, where the coupling terms appear. 

Proof. In order to prove identity (3.9) note that 

1 - 1 (9I(/)'^P 2 

-div.pf = -^Xm(0^A0f), = ^UemS^ 

Since 0? solves the corresponding equation in system (-FL), we can multiply the equation 
by 01 and add this identity to its conjugate; the conclusion follows from the properties 
of the nonlinearity. Concerning identity (3.10), observe first that, setting (pf)j(a;, t) = 
e^~^Xm{(t)^[x,t)dj(t)\{x,t)) for any j and dj = dx^, it holds 

= e'-''Xm{dt'^\d^<P\) + e'~''Xm{^,d,{dt<P\)) 
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In particular the second term integrates to zero. Concerning the first addendum, take the 
first equation of system (F^), conjugate it and multiply it by 2e~^dj(f)i. It follows 



2 



dx. 



p+1 J 'Kp 
Of course, one can argue in a similar fashion for the second component 02- Then, taking 
into account that all the terms in the previous identity but djV {x)]^)^^ and l^ll^'^^^^il^il'''''"^ 
integrate to zero due to the regularity of (pi, we reach 

,3.1) / Mk,. . I _ !| I 

Adding these identities for any j and taking into account that by the regularity properties 
of (pl it holds J dj{\(j)l\^~^^\(f)2\^~^^)dx = 0, formula (3.10) immediately follows. ■ 

4. Energy, mass and momentum estimates 

4.1. Energy estimates in the semiclassical regime. In order to obtain the desired 
asymptotic behavior stated in Theorems 1.5, 1.9, 2.2 and 2.4, we will first prove a key 
inequality concerning the functional S defined in (1.6). As pointed out in the introduction, 
the main ingredients involved are the conservations laws of the Schrodinger system and 
of the Hamiltonians functions and a modulational stability property for admissible ground 
states. 

The idea is to evaluate the functional S on the vector T'^ = {vf, i;|) whose components 
are given by 

(4.1) vl{x,t) =e-^«'W-["^+^iWl0f(ex + Xi(t),t) 

where X = (xi,X2), S = (^1,^2) are the solution of the system (H). More precisely, we will 
prove the following result. 

Theorem 4.1. Let $^ = (0f,0|) be a family of solutions of (F^), and let be the vector 
defined in (4.1). Then, there exist eo and such that for every s G (0,£:o) fl'^c? for every 
t e [0,T^), it holds 

(4.2) Q<E(T')-8{R) <a' + r]' + 0{e^), 
where we have set 

(4.3) a^{t) = \{Ut),Ut)) ■ {a\{t),al{t% r^%t) = W.it) + m\. 
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aj, rji are given in (3.5), (3.8) and R = (ri, r2) is the real ground state belonging to the class 
TZ taken as initial datum in {F^). Moreover, there exist families of functions 9f, y\ and a 
positive constant L such that 



(4.4) 



for every e G (0,5o) O'^d all t G [0,Tf). 

Proof. By a change of variable and Proposition 3.1, we get 



(4.5) 



Hi-Ml = \mex + x,it),t)\\l = — ||0f(.,t)||2 = m,, t > 0, 2 = 1,2, 



where nii are defined in (1.17). Hence the mass of vf is conserved during the evolution. 
Moreover, by a change of variable, and recalling definition (1.16) we have 



|V$^||2 + -(mi|eir + m2|6|') 



- / ■ {pl{x,t),pl{x,t))dx. 



Then, taking into account the form of the total energy functional, we obtain 



£(T") = E%t) 



[v{xm\^ + w{xm\'] dx+- (mii^ip + msi^r) 



Moreover, using Proposition 3.1 and performing a change of variable we get 



E'{t) = E'{0) = E'(ri 



X — X 



X — X 



^(i?) + -(mi|eir + m2|6r) 



+ J [V{ex + x)\ri\'^ + W{ex + x)\r2\'^] dx, 
this joint with Lemma 3.3 and the conservation of the Hamiltonians Hi{t) yield 



Mt),Ut)) ■ {Pl{x^t),pl{x,t))dx 
+ miV{x) + m2W{x) j [yi^Wi? + W {x)\(t)l\'^]dx 

- j {Ut).Ut)) ■ {pl{x,t),pUx,t))dx 

j [V{x)\<P\\' + W{x)\<PVi']dx + 0{e') 
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Using the definitions of and rji, we get 

- ^ y"(l - x{x)) [V{x)m\' + W{x)\<j>l\']dx + 0{e'), 

Since V and W are nonnegative functions, by (4.3) it follows that 

S{T%t)) - S{R) < a%t) + 7]%t) + 0{e^). 

Finally, (1.5) and (4.5) imply the first conclusion of Theorem 4.1, where the positive time 
is built up as follows. Let Tq > (to be chosen later). In order to conclude the proof of the 
result, notice that ai{t) and r]i{t) are continuous functions by Lemmas 3.5 and 3.7. Moreover, 
let rxe(f) be the positive number given in (1.9) for $ = T^(t). Notice that {t ^ ^r^{t)} is 
continuous and, in view of the choice of the initial data (1.3), it holds rxe(o) = 0. Hence, 
for every fixed ho, hi > 0, we can define the time > by 

(4.6) := sup{t G [0,T°] : max{a^(s), r]%s)} < ho, Tr^i^s) < hi, for all s G (0,t)}, 

Notice that, by (4.2) and choosing sufficiently small we derive, for all t G [0,T^) and 
e G (0,£o), that < £^(T^(t)) — £{R) < 3ho. Now we choose hi so small that hi < A, 
where A is the constant appearing in the statement of the admissible ground state (see 
Definition (1.3)). Therefore, from conclusion (1.10), there exists a positive constant L such 
that 

(4.7) Tr^^t)<Li£iT'it))-£{R))<L[a'{t)+r]'{t) + Oie')], 

for every t G [0,T^) and all e G (0,5o)- In turn, there exist two families of functions y'^(t) 
and 9l(t) i = 1,2 such that 

(4.8) ||T^(-,t) - (ei^'!Wri(. + r(t)),e^^'iWr2(- + r(t)))lle < L[a%t) + v%t) + Oie')] 

for every t G [0, T^) and all e G (0, Eq). Making a change of variable and using the notation 

et{t):=e9m, ym:=xi{t)-ef{t), 
the assertion follows. ■ 

Remark 4.2. The previous result holds for every t G [0,T^^) where is found in (4.6) and 
< Tq. But, we have not fixed Tq yet. This will be done in Lemma 4.6. 

4.2. Mass and total momentum estimates. The next lemmas will be used to prove the 
desired asymptotic behavior. We start with the study of the asymptotic behavior of the 
mass densities and the total momentum density. From now on we shall set 

a^{t) ■= a'{t) + \al{t) + a^(t)|, t > 0. 
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Lemma 4.3. There exists a positive constant Li such that 

/e^dx, / dx) -{mi,m2)6yn^t)\\{c^xc^Y 

+ \\P^{t,x)dx - M{t)5yi^t)\\(c^Y 

<L, [a-{t)+T^^{t) + 0{e^)], 
for every t G [0, T^] and e G (0, Sq). 

Remark 4.4. This result will immediately imply Theorem 1.5, once we have shown the 
desired asymptotic behavior of a^it) + r]^{t) and of the functional ^^e. 

Proof. For a given v G H^, a direct calculation yields 



l^'l j=l 



where Vj = v^^ and, in the last term, it appears the square of the modulus of the vector 
whose components are 1m{vvj). Then, we obtain 



In turn, using Theorem 4.1, it follows that, as e vanishes, 

\Im{vlVvl)\^ , f\Im{vlVvl)\^ 



,e|2 



0<£(N,h1)-£(i^) + / '^-(g^^^ dx^j 
<a'{t)+r^'{t) + 0{e^). 
Moreover, since |f||)||2 = || (^'f > ^^1) lb = II-RII2, we can conclude that 



dx 



(4.9) / J \ \ ^'^ dx + / J \ \ dx < a\t) + rt{t) + Oie"), 



for every e G (0,eo) and for every t G [0,T^^]. Using (4.1) and (4.5), for any i = 1,2 we get 

|Xm(^fV^f)P _ \eZm{^{ex + Xi, t)V0f (ex + Xi, t)) - (ga: + Xi, t)p|^ 
~ |0f(£a; + xi,t)|2 

2|Jm(0f(ex + Xi,t)V0f(ex + Xi(t),t))|^ 
|0f(ex + xi,t)|^ 
- 2EiiZm{(\)\{Ex + xi, t)V0- (ex + xi(t), t)). 

Whence, performing a change of variable and using definition (1.16), we derive 



I "^i I J I I 
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2 

dx + rrii 




'-'■1 









,£|2 



2 



e^mf 



/ pfda; 



which, by (4.5) is equal to (4.10). In turn, (4.9) imphes that 



(4.11) 



+ rrii 



f p^dx 



rrii 



_jv/2Pf(M) _ Ipfdx 10; 
10^1 mi 

<a'it)+r]'{t) + 0{e^). 

In order to prove the assertion, we need to estimate pf (t) for i = 1,2, where 



(4.12) pm 



4? / i^{x)\(f)'i\'^dx - miij{yl) 



P%x,t)ij{x)dx - M{t)i){y{) 



for every function -0 in such that 



C2 



< 1. From the definition of (3.5) it holds 



< 



< 



< 



P'{x,t)i){x)dx - M{t)i){y{) 
P'^{x,t)[ij{x)-ij{yl)]dx 
P'{x,t)[^{x)-ij{yl)]dx 



< 



i'iyl) ( I P%x,t)dx~M{t)^ 



+ \al{t) +al{t)\ 



^ — ^ rrii 



i=l 
2 



Pi{x, t)dx 



^(a:)|0f(a:,t)p 



dx - rriiiplyl) 
|0f(x,t)p 



dx 



for every e G (0,£:o) and for every t G [0,T^]. Taking into account that / Pidx is uniformly 
bounded and that, of course, 

-(yK(x,t)dx)^^_ 

there exists a positive constant Cq such that, if we set ip{x) = ■0(x) — ip{yl), it holds 

2 



pf(x,t) 



dx = 0, 



K^)<;^ / mx)\mx,t)\'dx+j2^ I \^{xm{x,t)\'dx 

|0f(x,t)P 



i=l 



E / 1^(^)1 



pl{x,t) 



rrii 



pI{x, t)dx 



-N 



dx 



+ \a\{t) + am\ + 0{e'). 
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From Young inequality and (4.11) it follows (from now on Cq will denote a constant that 
can vary from line to line) 

2 



+ 9E 



Co|^(x)| + i|^(a;)p 



\(t)l{x,t)\^dx 



1=1 



pf(x,t)£^/2 



|</)^(x,t)| rui 



pl{x,t)dx 



mx,t)\ 



+ \al{t) + am\ + 0{e') 
2 

i , 

< 



-y 



i=l 



1 



+ \alit) + al{t)\ + - [a^'(t) + r{t) + 0(e^)] , 
Using the elementary inequality a? < 26^ + 2{a — hY with 



1 /X — 
-r. ' 



and recalling that ip is a. uniformly bounded function we derive 



i=l 



Pt(i)<7NH / k|^(x)| + |^(x)| 



rff^Jll'Ux 



2 



+ \al{t) + ami + 2 k W + ^'(^) + ^(^')] • 

for every e G {0,eo) and for every t G [0,T^"]. Notice that satisfies the hypothesis of 
Lemma 3.3 and ip{yl) = 0, then by virtue of inequality (4.4) we obtain the conclusion. 



4.3. Location estimates for yl. In the next results we start the study of the asymptotic 
behavior of yf. 

Lemma 4.5. Let us define the function 



(4.13) 



I'it) = |7^(t)| + |72^(t)| , with Y,{t) = m^Xiit) - ^ I xx{xM{x,t)\^dx 



where x(x) is the characteristic function defined in (3.7). Then 'yl{t) is a continuous func- 
tion with respect to t and |7f(0)| = 0{e^) fori = 1,2. 

Proof. The continuity of 7"^ immediately follows from the properties of the functions x and 
01 . In order to complete the proof, note that Lemma 3.3 implies 



|7K0)| 



niiX 



[x + ey)x{x + ey)\r, 



'dy 



< Coe"^ + ImiX - xx{x)\\ri\\l2\ 
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and as x{^) = 1 we reach the conclusion. ■ 

Lemma 4.6. Let be the time introduced in (4.6). There exist positive constants ho and 
Eq such that, if \ri^ \ < ho and e G {0,eo) there is a positive constant L2 such that 

\xi{t) - ym < L, [a%t) + v%t) + 7^(t) + 0{e')] 

for every t E [0, T^] . 

Proof. First we show that there exist Tq > and B > such that 

(4.14) \ym\<B, 

for every t G [0,T^] with < Tq. Let us first prove that 

\\^yUt2) - ^vnh)\\c^* < B, for all ti, e [0,T^]- 

Let ^jj eC^ with ||^||c2 < 1 and pick ti, ta e [0,T^] with > ^i- From identity (3.9) and 
integrating by parts, we obtain 

^ J ^lj{x)m{xM)?-m{xM)?)dx = j^jyt j V;(x)M|l^cia; 

dt / il){x)d,i'v pl{x,t)dx 
< llV^lloo f'dtf \pt{x,t)\dx. 



It is readily seen from the estimate of V0f that the last integral on the right hand side 
is uniformly bounded, so that there exists a positive constant Co such that 

\\mx,t2)\Ve''dx-mx,t,)\ye''dx\\c^, <Co\t2-ti\<CiTo, 

with Ci = 2Co. Then Lemma 4.3, 3.5, 3.7 and 4.5 imply that the following inequality holds 
for sufficiently small e and ho (the quantity should be replaced by in the definition of 

^iW^lih) - ^l{h)\\c^' < CiTo + L[a'{t2) + r]'{t2) + a'{t^) + r^'it^) + 0{e')] 

<Cr[To + 0{e^) + ho]. 

Here we fix Tq and then So-, ho so small that Ci [Tq + 0{e^) + ho] < MKo where Ko is the 
constant fixed in Lemma 3.4 and from this lemma it follows 

\ylit2)-yl{h)\<C2Ko, 

for every ti, ta < ^o? and since ?/f(0) = x we obtain (4.14) for B = C2K0 + \x\. In view of 
property (4.14) we can now prove the assertion. Let us first observe that the properties of 
the function x imply 

\xiit) - yl{t)\ = — \mixi{t) - miyl{t)\ 



mi 

<—\im\ + — 

mi mi 



■^xxix)\(f)l{x,t)\^ - miyl{t) 
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Using (4.14) and (3.7) we obtain that xiul) = 1, so that there exists a positive constant C° 
such that 

\xi{t)-yl{t)\ < Co\\xx\\c4\<Pl\Ve''dx-m,Sy.\\c^*+Co'j%t). 
This and Lemma 4.3 give the conclusion. ■ 

In the previous Lemma we have fixed Tq such that also Lemma 4.3 and Theorem 4.1 hold 
and now we are able to prove Theorem 2.2. 

Proof of Theorem 2.2. We start the proof from the second conclusion of Theorem 4.1. 
By Theorem 4.1, the family of continuous functions g"^ : M"*" — M, 

(4.15) g%t)=L[a%t)+r]%t)+-f%t)] L = max{L, Li, L2} 
is such that g'^{0) = 0(£^) and it satisfies 

(4.16) 11$^ - (^ei(^«^+^^)ri(^^^^ , ei(^«^+^i)r2(^^^^ ^ £ < g%t). 
Moreover, Lemma 4.6 implies \ey^\ = \xi — yfl < g^{t), so that = Also, 



e 

2 



\r^{■) - r.(- - < ir I'll Vr.ll^. < C ( ^ ) , for alH = 1, 2. 



Then 

where we have set 

~g'{t) = g'{t) + C{g^{t)/ef. 
Since ^"^(0) = O^e"^), the assertion follows. ■ 

Proof of Theorem 2.4. In view of definition (4.15), the assertion immediately follows by 
combining Lemmas 4.3, 4.6 and 3.4. ■ 

4.4. Smallness estimates for d^, rj^, 7^. In the next lemma, under the assumptions of 
Theorem 1.5, we complete the study of the asymptotic behaviour of system (S'^) by obtaining 
the vanishing rate of the functions d^, r/^ and 7^ as e vanishes. The time To is the one chosen 
in the proof of Lemma 4.6. 

Lemma 4.7. Consider the framework of Theorem 1.5, that is V = W and C,i = ^2 = ^■ 

Then there exists a positive constant L such that 

a%t) +v%t) +7^(t) < L{T^)e^, for every t e [0,To]. 

Proof. By the definition of a'^{t) (see formula (4.3)) and taking into account that under 
the assumptions of Theorem 1.5 it holds C,i = C,2 = ^ (with respect to the notations of 
Theorem 4.1), there exists a positive constant C such that, for t > 0, 

a^t) = a%t) + \al{t) + a^(t)| < (1 + \m\)H{t) + a|(t)| < C|af (t) + a|(t)|. 
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Hence, without loss of generality, we can replace in the previous theorems (in particular 
Theorem 1.5) the quantities a^{t) and a'^{t) with the absolute value |af(t) + a2(^)|- In a 
similar fashion, it is possible to replace the quantity 7^(t) defined in formula (4.13) with the 
value |7f(t) + 72(^)1- We will prove the desired assertion via Gronwall Lemma, so that we 
will first show that there exists a positive constant L such that, for all t G [0,T^^], 



(4.17) 
(4.18) 
(4.19) 



&%t) < 0{e^)+L [\a%t)+7]'it)+^%t)]dt, 
Jo 

Jo 

I'it) < 0{e^) + L [ [a'it) + r]'{t) + ^'{t)]dt. 



Now, identity (3.10) of Lemma 3.8 yield 



d 
di 



«1 + «2)W 



< \\^V\\c4\<l>l\ -mi5^{t)\\c^' 

+ l|VV||c72|ll02lV^'^-"^25xwllc2* 



Hence, using Lemmas 4.3 and 4.6 one obtains, for all t G [0,T^^] 



d_ 
di 



< L^[a' + r]' + -f' + 0{e% 



for some positive constant Ai, yielding inequality (4.17). As far as concern r^^, using (3.9) 
and Lemmas 4.3 and 4.6 one has for t G [0, T^] that there exists a positive constant A2 such 
that, for all t G [0,T^ 



< 



mi VV(x(t)) ■ e(t) + m2Vy(x(t)) ■ ^(t) 
x{x)V{x)dw^pl{x,t) + / x{x)V{x)dw^p''2{x,t) 



[VixV)ipl+pl){x,t) - WVixit)) ■ {m,^{t) + m2m)]dx 
< \\V{xV) \\c^ \\P'{x, t)dx - M{t)6,^t) lie- 
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Let us now come to 7^. By the properties of the function Xi identity (3.9), Lemmas 4.3 and 
4.6 it follows that there exists a positive constant ^3 such that, for all t G [0,T^^ 



d f 

= j [V(a;x) ■ P'{x,t) - V{xx) ■ M(t)4(t)]rfa; 

< \\V{xx)\\c'^ \\P'{x, t)dx - M{t)5.,^t) ||c2* 

< A^[a' + r^' + + 0{e^)] 

Then inequalities (4.17), (4.18), (4.19) immediately follow from Lemmas 3.5, 3.7 and 4.5. 
The conclusion on [0,T^^] is now a simple consequence of the Gronwall Lemma over [0,T^^]. 
By the definition of and the continuity of a*^, a'^ and 77^ we have that = Tq provided 
e is chosen sufficiently small. To have this, one also has to take into account that, by 
construction (cf. formula (4.7)) and by the uniform smallness inequalities that we have just 
obtained over [0,T^^], we reach 



< L[a'{t) + ri'{t) + 0{e^)] < 0{e''), for all t e [0, T^] 



5. Proofs of the main results 

Proof of Theorem 1.5. In light of Lemma 4.7 we have f?^(t), ^^(t) < Le'^ for any t G [0, Tq]. 
Hence, the conclusions hold in [0, Tq] as a direct consequence of Theorem 2.2. Finally, taking 
as new initial data 

and taking as a new a guiding Hamiltonian system 

x{t) = e(t) 

|(t) = -VV{x{t)) 
x(0) = x(To), e(0) = e(To), 

the assertion is valid over [To,2To]. Reiterating (To only depends on the problem) the 
argument yields the assertion locally uniformly in time. ■ 

Proof of Theorem 1.9. Combining definition (4.15) with the assertions of Lemmas 4.3 
and 4.7, we obtain the property over the interval [0, Tq]. Then we can argue as in the proof 
of Theorem 1.5 to achieve the conclusion locally uniformly in time. ■ 
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